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Abstract 

Single molecule tracer diffusion studies of evaporating (thinning) ultrathin tetrakis-2-ethyl- 
hexoxysilane (TEHOS) films on silicon with 100 nm thermal oxide reveal a considerable slow¬ 
down of the molecular mobility within less than 4 nm above the substrate (corresponding to a 
few molecular TEHOS layers). This is related to restricted mobility and structure formation 
of the liquid in this region, in agreement with information obtained from a long-time ellip- 
sometric study of thinning TEHOS films on silicon substrates with 100 nm thermal or 2 nm 
native oxide. Both show evidence for the formation of up to four layers. Additionally, on 
thermal oxide, a lateral flow of the liquid is observed, while the film on the native oxide forms 
an almost flat surface and shows negligible flow. Thus, on the 2 nm native oxide the liquid 
mobility is even more restricted in close vicinity to the substrate as compared to the 100 nm 
thermal oxide. In addition, we found a significantly smaller initial film thickness in case of 
the native oxide under similar dipcoating conditions. We ascribe these differences to van der 
Waals interactions with the underlying silicon in case of the native oxide, whereas the thermal 

oxide suffices to shield those interactions. 
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Introduction 


Thin liquid films are widely present in nature, in particular as surface water layers. Also various 
technical applications make use of thin liquid films for lubrication, coatings or intermediate layers. 
Wettability considerations play an important role in the design of porous materials for catalysis and 
liquid storage. ^ Slippage, diffusion and flow characteristics of confined liquids have to be taken 
into account in nanoscopic Lab-on-a-Chip fabrication.^ Enhanced flexibility as well as economic 
considerations have led to an increasing interest in printing almost any kind of electronic devices.^ 
This also includes inkjet printing,^ for which droplet formation and evaporation is an essential 
feature. From numerous investigations, the surface roughness and chemistry is known to influence 
droplet shapes, wetting characteristics^ and thus also droplet evaporation.^ 

Deviations from the macroscopic behavior are observed on a molecular scale. For example, liq¬ 
uids show adsorption layers with a thickness of one or two molecular diameters. ^ Partially wetting 
liquids form a precursor film of similar thickness. ^ Seemann et al. found an influence of long-range 
capillary forces for nanometric polymer melt droplets in case of small contact angles (0 fa 5°). For 
larger contact angles (6 > 45°), the shape of the droplets could be described by surface properties 
and deviations invoking intermolecular interactions were only seen in the vicinity of the three- 
phase contact line.^ Investigations of thin liquid films by ellipsometry and X-ray reflectometry 
revealed molecular layering of up to four liquid layers at interfaces with solids, which was 
further confirmed by simulations. 

The question, whether the observed liquid layering has an influence on the molecular dynamics 
of the liquid, motivated optically detected single molecule tracer experiments in confined liquids. 
Thereby, a slowdown of the tracer diffusion in wetting precursor layers of liquid droplets, in ul- 
trathin liquid films and in a liquid film at a step edge on mica was observed. Furthermore, 
single molecule tracking (SMT) experiments, revealed a broad distribution of diffusion coeffi¬ 
cients, which indicates heterogeneous diffusion. The observation of heterogeneous diffusion 
was attributed to liquid layering. However, it has to be mentioned, that direct evidence of 

layering was obtained from investigation (by X-ray and ellipsometry) only for films on silicon 
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substrates with a thin (native) oxide layer. 

On the other hand, for sensitivity reasons, tracer diffusion studies were carried out only on 
silicon substrates with 100 nm (thermal) oxide, quartz glass or mica. Nevertheless, recent long 
term single molecule experiments on thinning liquid films on 100 nm oxide further addressed the 
mutable influence of molecular layering on liquid dynamics. They showed two discernible 
diffusion coefficients questioning the extend of the molecular layering to four layers as reported 
form ellipsometry and X-ray studies. 

These results initiated considerations about the influence of subsurface substrate properties on 
liquid dynamics and layering. Until now, no studies have been reported regarding the influence 
of the interface composition both on liquid layering and on diffusion dynamics close to the solid- 
liquid interface. However, for thin polymer films the wetting conditions are known to be influenced 
by subsurface material properties. Recent studies also demonstrate that the adhesion of proteins, 
bacteria and setal arrays of geckos on silicon wafers is influenced by the thickness of the silicon 
oxide interface. Additionally, in a previous study, we reported the effect of the subsurface 
material on the structure and dynamics of thin smectic liquid crystal films. 

The above mentioned tracer diffusion experiments on thinning TEHOS (tetrakis-2-ethyl-hexoxy- 
silane) films on silicon with 100 nm thermal oxide, also point to an influence of the subsurface 
substrate composition. However, for this type of multi-component interface no layering studies 
(neither by ellipsometry nor by X-ray) are known, which prevents from a direct comparison of 
diffusion and layering experiments. On the other hand, single molecule fluorescence experiments 
cannot be used to study diffusion in ultrahin liquid films on thin native oxide, since the small 
distance of the fluorescent tracer molecules to the underlying silicon leads to quenching of the 
fluorescence by electromagnetic interactions with silicon. 

To circumvent these experimental limitations, we will report optical ellipsometry experiments 
on thinning TEHOS films prepared both on 100 nm thick thermal oxide and on native oxide. 
Comparing the ellipsometry study with tracer diffusion experiments on thinning TEHOS films 
on thermal oxide, we describe the results in the frame work of long-range van der Waals 
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interactions. We find an agreement between reduced tracer mobility with decreasing thinning 
rates for decreasing film thickness and conclude that the multi-component interface properties also 
influence diffusion dynamics on a molecular scale. 


Experimental 

Substrate characterization and sample preparation 

(lOO)-polished n-doped (phosphor, conductivity 1...20 Qcm) silicon wafers with native oxide 
(Center for Microtechnologies, Chemnitz) and with 100 nm thermal oxide (CrysTec, Berlin) were 
rinsed alternately with acetone and ethanol (both technical grade) and sonicated for 10 min in ace¬ 
tone and afterwards in ethanol at 50 °C (both Merck, spectroscopic grade). In between and at the 
end the substrates were rinsed with de-ionized water and dried by nitrogen. AFM measurements 
showed smooth surfaces, with slightly larger average roughness for the thermal oxide 2.2 ±0.2 A 
compared to the native oxide 1.6 ±0.1 A. Macroscopic contact angles of water, diiodmethane and 
tetrakis-2-ethylhexoxysilane (TEHOS) (see 1) were similar for both types of oxide. Thus, the 
contact angles depend only on superficial properties of the substrates. 

Table 1: Macroscopic contact angles d [°] on Si wafers with different oxide thickness 


liquid 

0 on native oxide 

0 on 100 nm oxide 

water 

63±2 

61.4±1.3 

diiodmethane 

54.7 ±1.7 

52±3 

TEHOS 

5.6±0.9 

7.2±0.9 


The substrates were stored for 2 h under dry nitrogen atomsphere before preparation of the 
TEHOS films. Then, the substrates were dipcoated (KSV instruments) in solutions of 0.3% 
TEHOS (ABCR GmbH, Germany, for chemical structure see 1 c) in hexane (Merck, Germany, 
spectroscopic grade) using fast immersion, 20 s waiting time and 5 mm/min withdrawal speed. 
TEHOS is a branched flexible spherical molecule with a diameter of about 10 A. The chemical 
structure is given in 1 (c). 
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For diffusion measurements, Rhodamine B (RhB), Rhodamine 6G (R6G) and Oregon Green 
(OG), were alternatively used as fluorescent tracers, see 2 (b-d) for the chemical structures. The 
tracer molecules were added to the solution before dipcoating, resulting in nanomolar concentrated 
solutions in the ultrathin films. 

Diffusion measurements 

Single molecule tracking (SMT) experiments were carried out with a home-built wide-field mi¬ 
croscope described perviously, employing the 514 nm line of an Ar-ion laser, and an 

acquisition frame rate of 50 Hz, yielding a temporal resolution of 20 ms. 

For fluorescence correlations spectroscopy (FCS) a previously described home-built scanning 
confocal microscope was used with a pulsed excitation (40 MHz) at 465 nm, together with a 
hardware correlator (ALV-5000). This allows for detection of correlation times in the range of 
0.1 IJ.S to 10^ s. FCS curves were obtained from fluorescence intensity time traces lasting 300 s 
each. 

The film thickness for samples used with diffusion measurements was determined by varying 
angle ellipsometry with a VASE™ ellipsometer (J.A. Woollam Co., Inc., USA). The diameter of 
the ellipsometer spot was 3 mm. 

Mapping ellipsometry 

Mapping varying angle ellipsometry was conducted with a M-2000 ellipsometer (J.A. Woollam 
Co., Inc., USA), and spot diameter I mm. The samples were scanned on a rectangular grid with 
1 mm spacing, yielding 81 measuring points on a 10 mm x 10 mm sample. In a first step, the 
pure substrates were scanned. The obtained ellipsometry data sets were then used during further 
fitting to retrieve the film thickness of the prepared films. This step is necessary to correct for 
slightly varying silicon oxide thickness over the sample. The TEHOS film was modeled with a 
Cauchy model (Cauchy.mat, J.A. Woollam Co., Inc., USA) with the optical parameters A = 1.43 
and B = 0.002, which had been obtained before from a thick TEHOS film with the same instrument 
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(using the first two terms of the Cauchy power series for the refraction index n{X) =A + BX~'^ + 
Apart from ellipsometric measurements, the samples were stored in a chamber under 
a continuous flow of dry nitrogen gas at room temperature 21 ± 1°C. 

To obtain film thinning rates, succeeding ellipsometry scans were run on each sample in inter¬ 
vals of about 24 hours. Between succeeding scans an individual thinning rate r* = (hi — h 2 )/(At) 
was calculated for each measured point. Thereby, hi and h 2 are the film thicknesses measured in 
succeeding scans, and At is the time lag between these scans. The thus obtained thinning rates 
were binned according to a method described by Forcada and Mate: Each individual rate r* was 
added to all bins of width 0.5 A covered by the thickness interval from /i 2 to . In the end, the sum 
in each bin was divided by the number of its entries, thus yielding an average thinning rate r(h). 
Given a sufficient amount of data, this method enables the visualization of small changes in the 
thickness dependent thinning rate r(h), which otherwise would be hidden due to the long intervals 
between succeeding measurements. 

Results and Discussion 

The structure and dynamics of liquids at interfaces with solids can be studied in case of ultrathin 
liquid films via various experimental methods. Thereby, X-ray reflectometry will reveal the struc¬ 
ture of the liquid, whereas long-time studies on thinning films using ellipsometry are suited for 
investigation of the film morphology and slow macroscopic dynamics. On a microscopic scale, 
the dynamics within such films can be investigated via tracer diffusion experiments. Thereby, 
highly diluted fluorescent tracers are used to probe their local environments. Single molecule 
tracking (SMT) experiments are suitable for detection of diffusion coefficients up to the order of 
1 From trajectory analysis of SMT, detailed information about average diffusion 

coefficients Araj along trajectories and changes in the nature of diffusion (for example, due to con¬ 
finement to compartments)^^ is obtained. With respect to fast diffusion, the temporal resolution 
of SMT is more efficiently used (resolving diffusion coefficients up to about 10 jum^/s), when 
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analyzing probability distributions of diffusivities (scaled square displacements of detected tracers 
between sueeeeding frames). i’7d8,26 Considerable higher temporal resolution ean be aehieved by 
fluoreseenee eorrelation spectroseopy However, FCS is less sensitive for diffusion 

slower than 1 jum^/s.^^ We expect tracer diffusion dynamics to occur with diffusion coefficients 
between 55 /tm^/s (bulk value in TEHOS under ambient eonditions) and <1 /tm^/s as recently 

reported. i6>20,29 

We start our report with a long-time study of fluorescent tracer diffusion as a function of thin¬ 
ning in ultrathin liquid films on silicon with 100 nm thermal oxide. Thereby, we employ SMT with 
trajeetory as well as diffusivity analysis, and we use FCS on 100 A thiek films to amend our results 
from the SMT experiment on thinning TEHOS films on thermal oxide. A comparative study of 
traeer diffusion on thin native oxide is not feasible, due to the in that ease small distanee of the 
fluorescent tracer molecules to the underlying silicon (in the order of a few nanometers) which 
will lead to quenching of the tracer fluorescence. We therefore use ellipsometry for comparing 
the morphology and macroscopie dynamies of thinning films on both kinds of substrates and relate 
results about film structure and dynamics from tracer diffusion to those from ellipsometry. 

Tracer diffusion in ultrathin TEHOS films 

Previous work on tracer diffusion in ultrathin TEHOS films on silica substrates with 100 nm oxide 
revealed a considerable slowdown in comparison to bulk measurements by more than one order in 
magnitude. 16,19,20 information on the diffusion in dependence on the distanee to the solid- 

liquid interface could be obtained from a long-time study on thinning TEHOS films on silicon 
with 100 nm thermal oxide. Here we report an extension of this study, especially by FCS 
experiments on 100 A thick films, and a discussion in parallel to the results from the ellipsometry 
study, thus amending the previously published interpretations. 

FCS and SMT not only complement each other on the expected range of tracer dynamics in 
ultrathin TEHOS films. They also differ in sensitivity to adsorption events. If tracers will adsorb to 
the substrates, these events may dominate Araj- To select mobile tracers, a lower limit can be set to 
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the lateral extension (covered area) of trajectories. We name such trajectories "mobile." For long 
trajectories the average Djraj nevertheless may be small, if the dynamics are dominated by slow dif¬ 
fusion and adsorption. Fast diffusing tracers are likely to leave the detection area of about 250 /tm^, 
which renders short trajectories. However, determination of Dtj-aj from very short trajectories con¬ 
tains large errors.^® Within diffusivity analysis of SMT, immobile tracers may be excluded by 
setting a lower limit to individual diffusion steps. FCS is not sensitive to completely immobile 
molecules, since it is based on fluorescence fluctuations. However, ad-/desorption events during 
acquisition time will contribute and will result in apparent slow diffusion. By carefully comparing 
and combining results from these different methods, we obtain information about tracer mobility 
within ultrathin TEHOS films, which will be explained in the following and discussed further in 
combination with results from ellipsometry. 

Single molecule tracking experiments (SMT) using the tracer dye Rhodamine B (RhB, for the 
chemical structure see 2 d) follow changes in diffusion during film thinning over 294 h from sample 
preparation, starting with a 111 A thick film and ending with a 32 A one. 1 (a,b) shows the ob¬ 
tained broad distributions of diffusion coefficients Araj calculated via mean square displacements 
along identified trajectories. While 1 (a) collects all trajectories, 1 (b) puts the emphasis on 
mobility by selecting those trajectories, which exceed a lateral area of 1.3 [im^. The broad distri¬ 
butions of Dtraj seen in 1 (a and b) point to heterogeneous tracer diffusion, in agreement with recent 
reports. 1649,20 qualitative overall feature is that diffusion becomes slower with decreasing h. 

In principle the number of detected trajectories should be constant over the time of the film 
thinning experiment. However, the number of detected (mobile and immobile) trajectories in 1 (a) 
stays almost constant for h < 50 A, but is considerable higher for larger h. As was stated above, 
SMT, and in particular trajectory analysis of SMT, is restricted with respect to fast diffusion. The 
loss of fast tracers will artificially increase the number of detected trajectories, because missed out 
excursions into fast diffusion will split trajectories into several apparently separate trajectories (see 
supplement for further explanation). 

o 

Therefore, we conducted FCS experiments on 100 A thick TEHOS films within 24 h after 
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Figure 1: (a,b) Distribution of diffusion coefficients Dtraj calculated from trajectories of RhB in a 
thinning TEHOS film on 100 nm oxide. Indicated are the number of trajectories, mean values of 
Dtraj ( 0 ), the time after preparation in h, and the film thickness in A: (a) Dtraj from all detected 
trajectories (including i m mobile ones), (b) Dtraj for trajectories exceeding a lateral area of 1.3 [irn^ 
("mobile" ones only) and the ratio of these to all trajectories, (a,b) modified from. (c) Chemical 

structure of TEHOS. 


preparation using the Rhodamine derivates Rhodamine 6G (R6G) and Oregon Green (OG). The 
molecular structure of both dyes is depicted in 2 (b,c). Typical autocorrelation curves for both 
tracer molecules are shown in 2 (a). The autocorrelation function G 2 (t) for two-dimensional 
diffusion did not fit to the experimental ECS curves. A better approximation is realized with 
a two-component function Gbi(T), see ?? in the supplement. The fits are also indicated in 2 (a), 
fitting parameters are given in ??, diffusion coefficients in 2. Both dyes showed similar results 
within error, whereby the two components Di pcs = 0.1 ± 0.05 /s and .D2 ,fcs = 6 ± 5 ;Um^/s 

(averaged values from both dyes) were obtained with equal amplitude (see 2). A comparison with 
the distributions of Dtraj for similar film thickness (h = 111 A and h = 90 A, top two in 1 a,b) 
shows only some rare trajectories with Dtraj in the range of D2 ,fcs- Thus, it can be concluded 
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Figure 2: (a) Normalized FCS autocorrelation curves and two-component fits (- - OG, • • • R6G) in 
a 10 nm thick TEHOS film on thermal oxide, (b-d) Chemical structure of fluorescent tracers (b) 
Rhodamine 6G (R6G), (c) Oregon Green (OG), and (d) Rhodamine B (RhB). 

that a considerable amount of fast diffusion was missed out by the trajectory analysis of the SMT 
experiment. 

Now let us turn again to the distributions of Otraj shown in 1 (b). These trajectories can be 
considered as mobile, since they cover an area exceeding 1.3 jum^. We have already discussed the 
variation of the total number of trajectories (1 a) in relation to h. This explanation is also valid 
for the "mobile" trajectories. The number of detected mobile trajectories is largest for h^SO A, 

o 

and decreases with increasing as well as with decreasing film thickness. For h >90 A, the low 
number of detected trajectories can be explained by the above discussed loss of fast tracers. With 
decreasing h, the contribution from fast diffusion decreases considerably, increasingly preventing 
tracers to cover the area threshold set for mobile trajectories, while the total number of trajectories 
remains constant (la). Keeping these explanations in mind, the ratio of mobile to all detected 
trajectories can be used to explain tracer mobility. This ratio is indicated in 1 (b) and will be used 
later on for comparison with results from ellipsometry. 

o 

Below the film thickness of 32 A (about 3 times the diameter of TEHOS molecules), no "mo¬ 
bile" trajectories could be observed, which corresponds to the absence of mobile tracers up to film 
thicknesses h<10 A. It has to be mentioned that residual absorption layers and precursor films of 
wetting droplets are of similar thickness. In the precursor film of a TEHOS droplet on silica, tracer 
diffusion coefficients of R6G were reported to be in the range of 0.1 to 1 /tm^/s,^^ and thus 2-3 
orders of magnitude smaller than the bulk value T)buik ~ 50 ;Um^/s. 
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Improvement of temporal resolution with respect to trajectory analysis of SMT is gained by an¬ 
alyzing probability distributions of diffusivities (Jdiff = r^! (4t), i.e. scaled square displacements of 
detected tracer molecules at a fixed time lag t) for the smallest possible T = Tframe (the inverse ac¬ 
quisition frame rate). 17.18,24,26,32 analysis shows that the tracer diffusion can be approximated 
by two diffusion coefficients Di jiff = 0.3 ±0.1 /tm^/s and D? diff = 4.5 ±0.4 /tm^/s, which do not 
notably change in magnitude during film thinning, ^^48 ggg supplementary information for details. 
The T)/,diff are shown in 2. T)i,diff is within experimental error comparable to T)i,fcs reported for 
Alexa 488 by Grabowski et al.^® f^i.FCS for R6G and OG is subject to adsorption events which 
leads to smaller values. The fast components Z)2,diff and D2 ,fcs are comparable within experimen¬ 
tal error. A quantitative comparison shows that the amplitude a 2 ,diff of the fast component is only 
about 0.3 (see 2), less than the amplitude of Z)2 ,fcs (<72,fcs = 0-5). This supports our previous 
argument that some fast diffusing tracers are missed by SMT, even when using diffusivity analy¬ 
sis. The similar and T)2 ,fcs. show the influence of the strong vertical confinement on the 

observed diffusion coefficients, which prevents even PCS from resolving bulk-like tracer diffusion 
(f^bulk ~ 55 live?/ s). D 2 rather is an effective diffusion coefficient depending on physical diffusion 
coefficients and transition probabilities between the heterogeneous regions. For more details about 
the influence of heterogeneity and conflnement on diffusion coefficients ^^^2 ggg supplementary 
information. 

By using the cationic Rhodamine derivate R6G and the anionic OG for PCS experiments, 
we addressed the question, to which extend the diffusion of tracers is influenced by electrostatic 
(Coulomb) interactions with the substrate. As stated above, the diffusion coefficients are similar 


Table 2: Diffusion coefficients Di 2 [jUrn^/s] obtained from SMT and PCS on 100 A thick TEHOS 
films on thermal oxide. 


method 

tracer 

a\ 

Dx 

a2 

Di 

SMT (diffusivities) 

RhB 

0.7 

0.3±0.1 

0.3 

4.5 ±0.4 

PCS 

R6G 

0.5 

0.06 ±0.05 

0.5 

4±3 

PCS 

OG 

0.5 

0.14 ±0.05 

0.5 

8±4 

PCS 20 

Alexa 488 


0.2 


13 


error: standard deviation a (SMT), 0 + uncertainty of focal width (PCS). 
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for both tracer molecules, see 2. 2 also gives values for Alexa 488 in ultrathin TEHOS films 
on quartz substrates. In eontrast to RhB, R6G and OG, Alexa 488, does not adsorb to the sub¬ 
strate, thus leading to somewhat larger values Z)i,fcs = 0.2 /rm^/s and D2,fcs = 13 /tm^/s.^® 
The similar diffusion coefficients obtained for the differently charged tracer molecules, indicate 
that Coulomb interaetions are in the present ease negligible for traeer diffusion in ultrathin liquid 
films. Similarly, Hahl et al. found that Coulomb interaetions play only a minor role in protein 
adsorption. 

To summarize, the presenee of at least two different diffusion eoefficients indieates an inho¬ 
mogeneous film morphology. According to diffusivity analysis of the film thinning experiment, 
this heterogeneity lasts down to a film thickness of about 30 A. The slow diffusion coefficients 
determined via PCS and SMT can be assigned to diffusion processes in a near-surfaee region. Tra¬ 
jectory analysis of SMT indicates two significant changes in tracer dynamics and thus in liquid 
mobility. Above /j 80 A diffusion becomes very fast, indicated by the loss of fast trajectories. 
Below ~ 30 A diffusion beeomes considerable slow seen by the laek of "mobile" trajeetories. 
Since the fast component from diffusivity analysis is still present for h = 32 A (corresponding to 
not much more than 3 molecular diameters of TEHOS), the near-surface region of reduced mo¬ 
bility is limited to 1-2 diameters of TEHOS (about 10-20 A). ^^38 suggests an only weak 
layering structure of TEHOS on thermal oxide of probably only two molecular layers. In contrast, 
up to four liquid layers have been reported for TEHOS on silicon with native oxide. To examine, 
whether these findings are related to the film morphology, we compare the diffusion experiments 
with the ellipsometry study of TEHOS films on substrates with different oxide thickness reported 
in the following. A combined discussion of changes in mobility seen from tracer diffusion with 
features seen from film thinning rates will be given in the section on film thinning rates along 4 
(b). 
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Ellipsometry study of film thinning and morphology 


The recent development of simultaneous multi-wavelength ellipsometry together with an improved 
lateral resolution, facilitates long-time studies of thinning TEHOS films using mapping ellipsom¬ 
etry. In this section, maps of film thickness for both types of substrate are shown and discussed. 
Additionally, thickness dependent film thinning rates are calculated from the mapping ellipsometry 


data. 



Figure 3: Maps (10x10 mm^) of TEHOS film height z as a function of time (in h) after preparation 
on (top) 100 nm thermal and (bottom) native oxide. Arrows denote the direction of withdrawal 
during dipcoating. Maps are shown for thinning films from left to right. A false color code for the 
height z is given on the right. 


3 shows maps of TEHOS film thicknesses on 100 nm thermal (top row) and native (bottom row) 
oxide for two samples prepared under similar conditions. It is noticeable that the initially obtained 
TEHOS film on native oxide is thinner than the one on thermal oxide despite identical preparation. 
This is a general observation within our experiments. 3 (a and e) show thickness maps from the 
first ellipsometry runs on these samples 18 h and 20 h, respectively, after preparation. On 100 nm 
thermal oxide, the largest local film thickness was about 140 A, whereas on the native oxide the 
largest value (about 80 A) was about a factor of 1.8 smaller. We will explain these findings in a 
further section in terms of long-range electromagnetic forces. 

Both TEHOS films exhibited a convex surface across the sample. The partly observed elliptical 
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shape of both films is related to the preparation, since the dipcoating direction (indicated by arrows 
in 3 a,e) coincides with the small axes of the ellipses. The film on thermal oxide initially shows 
an inhomogeneous surface, which smoothens within the first 70 h. After this smoothening, the 
radius of curvature of the film surface on thermal oxide remain s smaller than that of the film on 
native oxide. This can be seen from the extension of the color steps in the map for the thermal 
oxide 260 h after preparation in 3 (b) in comparison to that in the map for the native oxide 66 h 
after preparation (f), which show similar similar film thickness. Further evidence for the different 
curvature is given by linear hight profiles, as shown in ??. Seemann et al. observed an influence of 
the interface potential on the shape of polystyrene nanodroplets in case of small contact angles.^ 
Our observation points to an influence of the interface potential on the shape of ultrathin film 
surfaces as well. Considerations and calculations of interface potentials will be postponed to a 
further section. 

Moreover, the film on native oxide shows a height maximum, which remains in the same po¬ 
sition (x ~ 6 mm, y ~ 4 mm) during film thinning (see 3 bottom row, and ?? -i- discussion in the 
supplement). This is not the case for the film on thermal oxide. Here a shift of the maximum 
is observed. All recorded maximum values and positions are given in?? in the supplement. This 
behavior corresponds to a flow of the liquid. The absence of TEHOS flow on native oxide is in 
accordance with the observation reported by Forcada and Mate, who found negligible flow at the 
edge of a TEHOS film generated by etching part of the film. Further information on the structure 
of the films can be obtained by calculating thickness dependent film thinning rates, which will 
be discussed in the following. 

Film thinning rates 

TEHOS films on native oxide show similar characteristics as those reported by Forcada and Mate. ^ ^ 
In contrast, flow of the liquid is observed for TEHOS only on thermal oxide. We used two different 
samples of TEHOS Aims on thermal oxide for comparison with our sample of a TEHOS film on 
native oxide, to examine how this finding depends on preparation conditions. 4 shows thickness 
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film thickness h [A] film thickness h [A] 


Figure 4: Film thinning rates for a TEHOS film on native oxide (A) and two different films on 
100 nm thermal oxide (•, o). (a) Thinning rates for the full data sets obtained over time spans of 
some weeks, (b) Thinning rates for one of the films on thermal oxide obtained after (V) one week 
and (o) more than four weeks, together with tracer mobility: (□) ratio of numbers of mobile to 
all trajectories, (■) amplitudes a 2 of fast components from diffusivity analysis - error bars denote 
standard deviations, and ( 0 ) amplitude of fast component from FCS. (c) Data shown in (a) on a 
magnified thickness range, arrows indicate minima of rate on native oxide (A), (d) Thinning rate 
for one of the samples on thermal oxide shown in (a, o) together with thinning rates calculated 
from several subsets (V, differing in colors) of the full data set (o). The subsets contain increasing 
numbers of succeeding measurements after sample preparation, and thus cover increasing time 
spans At of the total time of investigation. 


related thinning rates for the two films on thermal oxide (•, o) and the film on native oxide (A). 
After a period of about two weeks, the ellipsometric measurements on two of the samples (•, A, 
corresponding maps are shown in 3 c,b) had to be stopped due to experimental constraints. Thus 
the smallest studied film thickness is 13 A on native and 19 A on thermal oxide. The other 
thermal oxide sample (o) has been studied over a period of more than four weeks. 

o 

In general, the thickness dependent thinning rates r (given in A/h) are in a similar range for 
films on native and on thermal oxide. In particular, a constant r ^ 0.4 — 0.5 A/h is seen in 4(a) 
for a thickness interval between 90 A and 70 A, followed by an almost linear decrease to r ~ 
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0.1-0.05 A/h for/j < 40 A. Below about 40 A, the thinning rates decrease more slowly. Thus, 
two critical film thicknesses at ~ 80 A and ~ 40 A emerge, very similar to those previously 
discussed in diffusion experiments. The similar general range of r for films on native and on 
thermal oxide also agrees with the observation of similar macroscopic contact angles of TEHOS 
on both kinds of oxide (see 1), which are related to the liquid-vapor surface tension by Young’s 
equation. The latter is related to the boiling point of the liquid and thus to evaporation. For 
the film on native oxide (A) the constant thinning rate r starts at the largest observed local film 
thickness. Whereas, for the two films on thermal oxide, larger film thicknesses with much higher 
thinning rates are observed. 

The different regions of mobility on thermal oxide represented by the changes in thinning 
rate r are also seen from tracer mobility, as can be seen in 4 (b). For 100 A the ratio of 
mobile to all trajectories (□) as well as the amplitude of the fast component from diffusivitiy 
analysis of SMT (■) are low. The value obtained by FCS ( 0 ) is also included, indicating the loss 
of fast diffusion by SMT. This observation from tracer diffusion coincides with faster thinning 
rates (•, o) at > 90 A. The fast diffusion in this region probably is bulk-like, however, due to 
the strong vertical confinement the diffusion coefficient cannot be resolved by tracer diffusion, see 
supplementary information. For ~ 80 A, both, the ratio of mobile to all trajectories (□) as well 

O 

as a 2 (■) are high. The ratio of mobile to all trajectories (□) decreases for < 80 A, while a 2 (■) 
stays constant within error until 50 A. As was explained previously, the selection for mobile 
trajectories by an area threshold is rather arbitrary and in particular will miss out some mobile 
tracers at small h. The broad distributions of Araj from all detected trajectories even for = 32 A 
(see 1 a), points to considerable heterogeneous tracer diffusion down to this film thickness. Within 
the range of 80 to 40 A, the thinning rates obtained after one week (V) and after more than four 
weeks (o) differ considerably. As we will see later when discussing 4 (d), this can be attributed to 
the influence of liquid flow on the thinning rates on thermal oxide, which hampers the evaluation 
of thickness dependent evaporation rates. Nevertheless, the thinning rates show still considerable 

o 

liquid mobility for h down to 40 A. This agrees with the mobility from diffusivity analysis of SMT 
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(■). The change in mobility seen from thinning rates and from diffusivity analysis between 50 and 
40 A, might point to the onset of layering, since the extend of layering was previously reported to 
be about 4 liquid molecule diameters. 

For better visualization of potentially layer dependent thinning rates r, the data of 4 (a) are 
shown on a magnified thickness range in 4 (c). The thinning rate of the TEHOS film on native 
oxide (A) shows two minima at 23 and 33 A (indicated by arrows) with a separation of about 

o o 

10 A (see 4 c). A further minimum is expected for 13 A, which could not be reached due 
to experimental constraints. Nevertheless, the observation of local minima on native oxide with 
a separation matching the liquid molecule diameter agrees with the observation of Forcada et al., 
who attributed it to liquid layering. 

In contrast, for the films on thermal oxide (4 (c) •, o) no clear local maxima and minima can 

o 

be seen. The thinning rate for one of the films (•), shows a slight undulation for h <30 A. But, 
for this sample, the end of measurement renders low statistics in the range below 30 A. The other 
sample on thermal oxide (o), shows no undulations of the thinning rate for < 40 A. For this 
sample, the film thickness does not decrease below 20 A, although the sample was studied over a 

o 

period of more than four weeks. For < 45 A, the thinning rate r{h) of this sample is much lower 
than those of the other two samples. Yet, for TEHOS on thermal oxide, lateral flow is observed, 
which is expected to influence the thickness dependent thinning rates. 

Let us assume flow of liquid from position A to B between succeeding measurements. Then 
the local thinning rate for A appears to be higher than the average value for the corresponding 
thickness in absence of flow, whereas the rate at B is smaller than the average rate corresponding 
to the thickness at B. A flow from thicker regions to lower regions thus will explain the low thinning 
rates observed for this sample. If liquid flow is responsible for the here observed low thinning rates, 
a change in thinning rates should be seen, when average thinning rates are calculated from different 
time intervals At of the more than four weeks total period of measurement. 

4 (d) shows such thinning rates calculated from several subsets (V, subsets differ by color) for 
the film on thermal oxide shown in 4 (a-c, o). Thereby, all subsets start from sample prepara- 
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tion and contain an increasing number of succeeding ellipsometry scans, thus, spanning increasing 
time intervals At of the total period of measurement. The arrow marks the increasing time span At 
covered by the subsets. The full data set is denoted by (o). Obviously, for the thickness interval be¬ 
tween 80 and 20 A the apparent film thinning rates decrease with increasing At. This observations 
can be explained by a flow of the liquid from thicker regions as well as laterally outside regions 
into the region of measurement, whieh is reasonable, sinee also the edges and the baek side of the 
sample are covered by TEHOS. 1150 h after sample preparation, still a residual TEHOS film of 
20 A thickness is observed. 

From our film thinning experiments under ambient eonditions in air, we generally observed 
residual layers of thicknesses between 10 to 20 A. The actual residual thickness depends on chang¬ 
ing ambient eonditions. For the film on native oxide, no changes over time were seen in the cal¬ 
culated thinning rates (see ?? (a) in the supplement). For the other thermal oxide (4 a,c, •), the 
change over time of the apparent thinning was less pronounced (?? (b) in the supplement). For this 
sample, liquid flow is observed also from evaluation of height profiles (see 3 e-d, and supplemen¬ 
tary information), which contributes to the apparent thinning rates. According to our observation, 
calculated thinning rates on thermal oxide contain contributions from evaporation and from liquid 
flow at any particular interval of the total measurement time. However, on this mieroseopie scale, 
the extend and direction of flow is induced by thermal fluctuations and defects. Thus the ratio of 
both components contributing to the thinning rate may change over time. On thermal oxide, the 
contribution from liquid flow to the apparent thinning rates smears out any dependenee on layer¬ 
ing, whereas on native oxide the absence of liquid flow ^ ^ enables the observation of discrete steps 
related to layer dependent thinning rates. This is further supported by the larger radius of curvature 
(see previous section) for the film on native oxide, where the almost flat film surface faeilitates the 
observation of changes in thinning rate due to layering. 

o 

For the film on native oxide, molecular layering can be seen for thicknesses below 40 A, in 
accordanee with the report by Foreada and Mate. However, from our experiments, it eould not 
be determined, to which extend molecular layering is present in TEHOS films on thermal oxide. 


18 



Nevertheless, our ellipsometry studies reveal a thinner initial film thickness on native oxide and 
a higher mobility of TEHOS on thermal oxide. The larger radius of curvature for the film on 
native oxide, is in agreement with observations by Seemann et al., who found an influence of 
long-range capillary forces on the shape of nanometric polymer melt droplets in the case of small 
contact angles (0 Since, both kinds of oxides displayed similar macroscopic contact angles 

(see 1) and surface roughnesses, the observed differences for the liquid Aims cannot be explained 
by superficial properties of the substrates. However, long-range van der Waals (vdW) interactions 
may account for the observation, because they are different for both kinds of substrates. Thus, we 
will now turn to evaluate the contribution from vdW forces to the interface potentials for the here 
studied ultrathin TEHOS films. 


Interface potentials 

Van der Waals forces arise from electromagnetic fluctuations coupling between different (or like) 
materials across a medium. The direction and strength of the forces result from the existence or 
absence of resonances of these fluctuations in the relevant materials. VdW interactions combine 
zero frequency fluctuations of permanent dipoles with pure dispersive fluctuations (of induced 
dipoles only), see schematic presentation in 5. These forces are non-additive, apart from diluted 
gases. In 1956, Lifshitz applied a continuum approach in quantum field theory to describe vdW 
forces.His theory was simplified using a summation of oscillator frequencies.^'^ However, its 
application acquires extended knowledge of optical material properties, which are hardly available 
for solids, and which are even harder to achieve for soft matter. Nevertheless, an approximation 
may be obtained using the dielectric properties of the relevant materials in the visible frequency 
range. 

A thin liquid film of thickness on a homogeneous substrate in air corresponds to a 3-layer 
system, (see 5, left), for which the effective vdW potential is given by^^“^^ 


*^vdW ih) 


^12-23 


( 1 ) 
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Figure 5: Scheme of van der Waals interactions for (left) a TEHOS film with thickness h on thick 
thermal oxide (3-layer system, since the underlying silicon is shielded by the thick oxide), and 
(right) a TEHOS film with thickness h on silicon with a thin native oxide of thickness d (4-layer 
system). 


where A 12-23 is the material dependent effective Hamaker coefficient (details in the supplement). 
In case of an additional coating material on the substrate, ?? has to be extended including the 
interaction of the additional layer with the liquid film. Eor the system of a thin liquid film of 
thickness h prepared on a silicon wafer with an oxide layer of thickness d (see 5, right), this results 
in 


^vdw{h) 


^SiOl _ Asi-Asi02 
llTth^ nn{h + dY' 


( 2 ) 


where A 5,02 and are the effective Hamaker coefficients for the system Si02/film/air and Si/film/air, 
respectively. This approximation yields an one oder of magnitude stronger effective interface 
potential for a TEHOS film on native oxide in comparison to a TEHOS film on thick thermal ox¬ 
ide. The difference is caused by the influence of the strong polarizability of the underlying silicon, 
which is shielded in case of thick thermal oxide, but contributes in case of only 2 nm native oxide 
(5). A more detailed description is given in the supplementary information. 6 (a) shows the thus 
obtained effective interface potentials <l>(/i) for TEHOS films on both kinds of substrate. 

VdW potentials calculated from this simplified Lifshitz theory and its approximations are de¬ 
rived from (homogeneous) continuum material properties and thus do not account for modifications 
of these properties on a molecular scale (for example due to molecular layering). In particular, the 
refractive index is related to material density. Therefore, density oscillations derived from X-ray 
reflectometry in thin liquid films close to solid-liquid interfaces can be used to improve the de- 
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Figure 6: Effective interface potentials 4>(/?) for TEHOS films on silicon wafers with (—) native 
and (• • •) 100 nm thermal oxide (a) calculated from continuum theory ??, (b) with oscillations from 
layering, (c,d) compared to film thinning rates on (c) native and (d) thermal oxide (shown in 4 a). 


termination of the effective interface potential. Yu et al. studied molecular layering of TEHOS 
films on silicon with native oxide. They found electron density oscillations with a period of 11 A, 
which corresponds approximately to the molecular diameter of TEHOS molecules. The ampli¬ 
tude of these oscillations decayed within several tens of nanometers. Thus, their findings resemble 
a liquid layering up to four molecular diameters above the substrate, whereby, the localization of 
the molecules within the layers decreases with distance from the substrate. Yu et al. showed that 
the experimental data could be modeled by a damped sine function. Eor this, they used a model 
independent fit as well as a damped sine function to fit the data. Eor a distance larger than 10 A, 
both gave very similar results. Only very close to the substrate, the electron density cannot be 
reproduced by the damped sine function. Thus, we use a sine function of the same period (11 A) 
to modulate the interface potentials obtained from continuum material properties. The resulting 
potentials are shown in 6 (b). Although X-ray measurements are reported only for TEHOS on 
native oxide, a similar modulation is also superposed on the effective interface potential for ther¬ 
mal oxide. However, it is not straightforward to decide about the exact position of the absolute 
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phase due to possible conformational changes of the branched and thus soft TEHOS molecules in 
the layer next to the substrate. Due to this lack of information about the thickness of the first 
molecular layer, the phase shifts of the sine functions are chosen in such a way that the minima of 
the thinning rates at fa 30 A (thermal) and 35 A (native) are met by maxima of the potentials, 
as can be seen in the comparisons in 6 (c,d). As was shown by Yu et ak, the damped since function 
cannot represent the interface potential below 10 A. Therefore, the modulated potentials are only 

o 

considered for the discussion for h> 10 A. 

The differing strength of the interface potentials for ultrathin TEHOS films on native and on 
thick thermal oxide could already be seen from continuum theory (6 a). Yet, from the comparison 
of the layer-modulated potentials (6 b-d), similarities and differences can be seen even more clearly. 
Both potentials are large for = 10 A, and thus support a residual (frozen) adsorption layer. The 

o o 

potential maxima at ~ 20 A and 30 A seen for thermal oxide are much weaker as the maxima 

o 

for native oxide in the range below = 40 A. This may explain the higher curvature on thermal 
oxide and the flattened surface of the ultrathin TEHOS fllm, as well as the restricted macroscopic 
fluidity on native oxide. 

o 

Despite the differences in potential strength, the thinning rates for = 40 A are similar on 
both substrates. We found that the thinning rate for thermal oxide is just an apparent thinning 
rate influenced by evaporation as well as by liquid flow, and thus the evaporation rate may even 
be higher, as is seen in 4 (b) for thinning rates obtained from shorter spans of measurement time 
At. Nevertheless, in the range between 40 and 70 A, the thinning rates increase to a constant 

o 

value r 0.4 — 0.5 A/h. This points to geometrical constraints influencing the film thinning rates, 
in particular, molecular layering may also extend to a four molecules distance on thermal oxide. 
This suggestion is further supported by the drop in the fast component of tracer diffusion observed 
at this film thickness (see 4 b, ■). Our calculations did not take the finite velocity of light into 
account, which leads to retardation of the interaction for a film thickness h> 50 A, implying an 
even faster decay of the interface potentials with The increase of the thinning rate r{h) 

starting at 40 A, even for native oxide, points to a general crossover from an evaporation rate 


22 



mainly influenced by $(/?), to an evaporation rate influenced only by the surface tension of TEHOS 
in air. 

The stronger potential in general, does not necessarily explain the observed smaller over¬ 
all film thickness from preparation on native oxide. Here, a closer look at the contributions to 
the potential may provide further information. 6 shows effective interface potentials containing 
contributions from van der Waals interactions related to both, permanent dipoles (zero frequency 
fluctuations) and dispersion, i.e. fluctuations from induced dipoles. The contribution of the fluctua¬ 
tions from permanent dipoles to the effective Hamaker coefficients is negative in case of the system 
Si/TEHOS/air (enhancing film stability), but it is positive for the system Si02/TEH0S/air (see ??in 
the supplementary information). A similar trend is observed for systems containing hexane instead 
of TEHOS (see supplement). 

Silicon wafers are cut from highly ordered crystals. Therefore, its permanent dipoles are also 
highly ordered. In technical crystal grow, the so called "van der Waals epitaxy" uses ordering prop¬ 
erties of van der Waals forces to grow crystals on substrate crystals in a highly ordered structure 
in spite of a mismatch in lattice properties. Atomistic simulations of liquids at interfaces with 
crystalline solids showed substantial in-plane ordering within liquid layers in case of compatible 
structures of the liquid phase and the solid. For incompatible liquid phases the liquid formed lay¬ 
ers with a liquid-like in plane structure. It is reasonable that not only solids, but also liquids are 
influenced by the crystal structure of the underlying substrate crystal in the range of several tens 
of nanometers, which is likely to lead to restricted rotational and translational motion of liquid 
molecules in this interface region. By this, the restricted mobility and highly ordered structure of 
ultrathin TEHOS films on Si wafers with native oxide can be explained. Furthermore, the strong 
ordering and reduced mobility may allow attachment of TEHOS molecules only at preferable sites 
within the growing film during dipcoating, and therefore may account for the thinner TEHOS films 
on native oxide compared to those on 100 nm thick thermal oxide. This suggestion is further sup¬ 
ported by the recent observation of the coordinated mobility of aggregated fullerenes on weakly 
interacting substrates. 
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The observed flow of TEHOS on the thick thermal oxide may be related to thermally induced 
thickness fluctuations and thus to capillary waves at the free surface with air. In case of thin poly¬ 
mer films, surface modes with growing amplitudes are responsible for spinodal dewetting. These 
growing modes will occur, if the interface potential contains a local minimum and a negative sec¬ 
ond derivative around a critical film thickness The interface potentials for the here 

studied TEHOS films, however, do not contain such minima (since sin‘‘{x) = —sin{x)). Yet, the 
surface of the TEHOS films on thick thermal oxide was considerably convex. Thus, the derivation 
given by Vrij for the increase of free energy AF associated with thickness fluctuations of a film of 
otherwise homogenous thickness does not fully apply. In particular, the contributions from the 
local thickness modulations to AF will override any contribution from the considerably weak 
for film thicknesses h>20 A. 

In hydrodynamics, liquid flow at interfaces of solids is discussed in terms of slip or no-slip 
boundary conditions. Thereby, the occurrence of slip is usually related to superficial properties 
like interface roughness and chemistry. The above reported results emphasize that also in this case 
long-range interactions with subsurface materials have to be considered. According to theoretical 
considerations by Bocquet and Barrat, a reduction of will enhance the slip length and thus 
the flow of the liquid, which is in accordance with our above described observations. The ob¬ 
served impact of the interface potential on liquid dynamics, may also contribute to the mobility 
of polymer chains at interfaces even besides dewetting situations. Thus recently observed devia¬ 
tions of polymer properties in thin films from those in the corresponding bulk materialshould be 
discussed also in the light of these long-range vdW interactions. 

Conclusions 

This contribution aims at the understanding of the influence of interface properties on formation 
of structure and molecular mobilities in liquid films on a nanometer scale. The observations from 
film thinning of ultrathin TEHOS films on silicon with either native or 100 nm thermal oxide 
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show different behavior of the films on the two kinds of substrate. In accordance with reports 
from literature, TEHOS films on native oxide are highly structured and show negligible flow. 
In contrast, on thick thermal oxide flow was observed for TEHOS films down to a residual film 

o 

thickness /j fa 20 A. This flow may conceal layering from detection by our ellipsometry study. 
Nevertheless, film thinning rates showed an overall similar thickness dependence on both kinds of 
substrates. In particular, the thinning rates r decreased from r fa 0.4 A/h at fa 70 A to r fa 0.1 A/h 
at fa 40 A. This points to a similar structural change on both kinds of substrates in the region of 
fa 40 A. Despite the higher mobility and surface curvature for films on thermal oxide, TEHOS 
may be organized into molecular layers also in this case. This suggestion is further supported by 
our findings from tracer diffusion in thinning TEHOS films on thermal oxide. These experiments 
show changes at fa 80 A (for thicker films, faster diffusion occurs), and fa 50 A (number 
of detected trajectories decreases to a constant value and amplitude a 2 of fast component starts 
decreasing with h). The latter corresponds to the expected range of molecular layering. Further 
information may be obtained from a comparative X-ray study. 

The here reported observation of subsurface composition influencing both the structure and dy¬ 
namics at solid-liquid interfaces, points to the necessity of taking also long-range vdW interactions 
into account. Pairwise additivity of Lennard Jones interactions may overestimate the extend of lay¬ 
ering reported from simple Monte Carlo simulations. In a similar way, for modeling adsorption 
of organic molecules on metals it was found necessary to include vdW interactions in DFT calcula¬ 
tions.^^ From our findings, it can be concluded that the structure and mobility of liquids in vicinity 
of solid-liquid interfaces is influenced not only by superficial properties of the substrate, but also 
by its subsurface composition. In particular, liquid structure formation and diffusion dynamics 
observed for one specific kind of substrate, may not be generalized to all substrates containing 
similar surface chemistry and roughness, while neglecting mutable influence of different subsur¬ 
face composition. On the other hand, this very influence of the subsurface material may be used 
for tailoring structure and dynamics within technical applications without the need of changing the 
superficial properties of the substrates. 
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Influence of confinement and heterogeneity on tracer diffusion coefficients 

In case of uitrathin films, neither single molecule tracking (SMT) nor fluorescence correlation 
spectroscopy (PCS) can resolve the vertical (z-)component of diffusion, but will detect only two- 
dimensional diffusion (the focal depth is about 1 jim and thus at least two orders of magnitude 
larger than the film thickness h). We will explain the influence of film heterogeneity and confine¬ 
ment on the observed diffusion coefficients using a simple two-layer system with layer-dependent 
diffusion coefficients D\f^D 2 and exchange rates /12 and / 21 , where /12 is the exchange rate from 
layer 2 into layer 1. Then, equilibrium probabilities pi to find the dye in layer i are given as^’^ 


Pi = 


fi 


ij 


fij + fji 


*To whom correspondence should be addressed 


(SI) 


1 



In the long time limit of the observation, the effective diffusion coefficient of the system can 


be calculated from the pi and the corresponding diffusion coefficients Di 


D^n = ^PiDi, (S2) 

i 

which is a weighted sum of the diffusion coefficients related to the involved heterogeneous re¬ 
gions. The observed diffusion coefficients from the two layer system depend on the relation of 
the exchange rates fn and /21 to the observation time Tobs-^~^ For Tobs ^ I //12 ^nd I// 21 , the 
experiment will yield In case Tobs is much smaller than both inverse exchange rates, the layer 
dependent diffusion coefficients D\ and D 2 will be observed with the respective amplitudes pi and 
P 2 . In between these two cases, the temporal behavior of the system is complex and cannot be 
predicted easily. 

The diffusion coefficient of the used Rhodamine dyes in bulk TEHOS is about 55 jum^/s.^’^ 

o 

Assuming a homogeneous TEHOS film with thickness h= 100 A and no structural change through¬ 
out the film, it is possible to evaluate the vertical mean first passage time (mfpt) for tracer molecules. 
For a reflecting boundary condition at the liquid-air boundary and an adsorbing boundary condi¬ 
tion at the substrate, the vertical mfpt is shorter than 10“^ s for Z) > 5 ;Um^/s.^ In other words, a 
tracer molecule desorbing from the substrate will take on average less than 10“^ s to be resorbed 
to the substrate under the assumption of a probability of 1 for resorption. ^ 

The temporal resolution of ECS is in principle suited to detect time as short as 0.1 ps = 10“^ s. 
Although, the obtained autocorrelation curves from the PCS experiment did contain components 
of that order, ^ the ECS curves were only analyzed down to 0.1 ms here. So the vertical mfpt 
could not be resolved directly. Nevertheless, the fast vertical passage of a 100 A thick TEHOS 
film has another influence on obtained diffusion coefficients. With PCS, translational diffusion 
coefficients are obtained from passages of the tracers through the focal area. The diameter of 
the focal area in the here used setup is 0.8 jum. Using the Einstein-Smoluchowski Equation D = 
(r^)/(4T) for two-dimensional diffusion,^ the mean square displacement during the mfpt can be 
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calculated as (r^) = ADx = 2.2 • 10“^ iim^. This is considerably smaller than the focal area. 
In a 100 A thick heterogeneous film, the diffusion coefficient obtained by PCS, for this reason, 
will be an average diffusion coefficient for the passage of the focal area. To resolve physical 
diffusion coefficients from the observed ones, the vertical structure of the film and the exchange 
rates between heterogeneous regions have to be known, which up to now is not the case. 

With SMT, the smallest possible observation time is the inverse frame rate Tframe> which in the 
here used experiment is 20 ms. This is several orders of magnitude larger than the above evaluated 
mfpt for a homogeneous film of thickness = 100 A. As a consequence also here the observed 
diffusion coefficients are strongly influenced by the vertical confinement. In particular, a bulk-like 
layer with diffusion coefficient D fa 50 ;ttm^/s cannot be resolved. Since Tobs S> mfpt even for a 
100 A thick layer, it will not make much difference on the Deff, whether the tracer can diffuse with 
Dbuik in a 20 A thick layer or in a 100 A thick layer. Due to the fast mfpt, the transition rate out 
of this layer will mainly depend on the probability to absorb into the region with smaller mobility, 
and not on the thickness of the bulk-like layer. 

Araj obtained from trajectory analysis of SMT will suffer further bias from heterogeneity. Fast 
tracer molecules will be lost to some extend during tracking. From the evaluation of trajectories 
we know that tracer molecules will experience periods of slow and of fast diffusion within ultrathin 
TEHOS films. When a tracer is lost from tracking during a fast period, but can be tracked before 
and after this excursion into fast diffusion, the SMT experiment will show two separate trajectories 
(with small diffusion coefficient) instead of a single longer one. This one the one hand will in¬ 
crease the number of detected trajectories. On the other hand this will overestimate slow diffusion, 
because there now will be two trajectories with underestimated diffusion coefficients contributing 
to the distributions of Dtraj, while the lost period of fast diffusion contributes to neither one of 
them. As long as only some short events of fast diffusion are lost, the remaining distribution of 
Araj will show a considerable amount of fast trajectories. If nearly all excursions of tracers into 
fast diffusion will be missed out, this will render only few remaining trajectories of considerable 
length which exceed the set lower limit on covered area to be taken into account for "mobile" Dtj-aj. 
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Diffusion coefficients derived by FCS 

Due to the small focal area, the exciting laser power is considerable high in case of FCS exper¬ 
iments. Within our experiment it was 0.4 kW/cm^ (compared to 0.1 kW/cm^ for SMT). This 
results in enhanced photobleaching of the dye molecules. In particular, RhB used within SMT 
experiments bleached fast, yielding FCS curves within noise. R6G and OG were more stable, 
rendering reasonable autocorrelation curves as can be seen in ??. 

For fitting the FCS autocorrelation curves for translational diffusion, a two-dimensional ge- 

o 

ometry is assumed, since the film thickness of 100 A is much smaller as the focal depth of the 
illuminating laser beam/detection area. A single component correlation function for diffusion in a 
two-dimensional geometry did note fit to the experimental data. Fitting could be improved by 
using the following two-component correlation function 




(S3) 


From the characteristic time Td related to translational diffusion the diffusion coefficient D can 
be derived by^’^ 


^xv 

D= ^ 


4to’ 


(S4) 


where Wxy is the lateral focus width (from center to 1/e decay of intensity) of the illumination/detection 
area on the sample, which was 0.4 ;Um in our case. Parameters from fits to the FCS curves using 
?? are shown in SI . 

Table SI : Fitting parameters of FCS autocorrelation curves according to ?? and derived diffusion 
coefficients using ?? for 100 A thick TEHOS films. 


tracer 

a\ 

Tz)i [ms] 

Di 

[S] 

D2 

R6G 

0.5 

9.4 

0.06 

0.28 

4 

OG 

0.5 

5.0 

0.14 

0.70 

8 
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Diffusion coefficients obtained from probability distributions of diffusivities 

Diffusion coefficients Dtraj calculated from mean square displacements along trajectories as shown 
in ?? are strongly averaged. In particular, fast diffusion will not be resolved by this method. 
Probability distributions of diffusivities Jjiff ( /{Ax), i.e. scaled square displacements 
during a time lag x) are superior in respect to this and they also yield larger data sets.^’^ 





Figure SI: (a-d) Probability distributions C(Jdiff)T: of diffusivities rfdiff from long term observation 
of RhB in a thinning TEHOS film, for film thickness h\ (a) 111 A, (b) 90 A, (c) 78 A and (d) 53 A. 
The 10 different data sets for each film thickness are plotted with different colors. Dashed lines 
show examples of bi-exponential fits, (e) Diffusion coefficients Di jiff, D 2 ,diff from bi-exponential 
fits to probability distributions of diffusivities for each data set for decreasing h from 111 A to 
32 A. Mean values are connected by lines, (f) Corresponding amplitude a 2 of the fast component 
^2,diff- Error bars denote standard deviations. (e,f) modified from.^’® 


In SI (a-d) probability distributions of diffusivities are shown corresponding to the film thin¬ 
ning experiment shown in ??, for film thicknesses (a) 111 A, (b) 90 A, (c) 78 A and (d) 53 A. For 
the remaining five investigated film thicknesses between 49 A and 32 A, the corresponding prob¬ 
ability distributions were previously published.^’® The C(^/diff)T contain only mobile ^/diff- Diffu¬ 
sivities related to completely immobile dye molecules had been excluded using a threshold related 
to the signal to noise ratio of the respective detected dye molecules.^’* This enables fitting the 
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C(Jdiff)T with a bi-exponential function according to^’^ 


C(i/diff)r = «iexp -Fa2exp ai+a2 = l. (S5) 

V ^l,diff/ V ^2,diff/ 

Examples of such fits to C(d?diff)T are indicated as dashed lines together with the derived diffusion 
coefficients in SI (a-d). SI (e) shows the obtained diffusion coefficients from fits to C(Jdiff)T for 
h decreasing from 111 A and 32 A using ??. Mean values of Z),,diff derived for each value of h 
are connected by lines. They are almost constant during film thinning, yielding a fast component 
f^ 2 ,diff = 4.5±0.4 ;Um^/s and a slow component T)i diff= 1.0±0.2 jum^/s. The latter is influenced 
by the for this analysis set threshold for mobile tracers, which was chosen to ensure exclusion 

o 

of all immobile tracers. Analysis of further SMT experiments on 60-80 A thick TEHOS films 
employing tri-exponential fits on probability distributions of diffusivities including immobile ones, 
yields Z)i,diff = 0.3 ±0.1 /tm^/s.^ In SI (f) the mean amplitudes a 2 are shown for decreasing h 

o o 

from 111 A and 32 A. Error bars denote standard deviations. The smaller a 2 for the two initial film 
thicknesses h = III A and h = 90 A can be explained by incomplete detection of fast diffusion, 
as is seen in comparison with results from ECS experiments on equal film thickness (see SI ). 

o 

The situation is different for h <50 A. There the smaller and almost constant number of detected 
trajectories (?? a) points to a complete detection of tracer molecules. The decreasing values of a 2 , 
therefore, are caused by a reduced mobility of tracer molecules, as is seen also by the change in 
the number of detected mobile trajectories in ?? (b) and in the ratio of mobile to all trajectories in 
?? (b). 

Liquid flow 

S2 gives the maximum thicknesses for the thinning TEHOS films on native and thermal oxide, 
for which maps were shown in ??. For the film on thermal oxide, the position of the maximum 
film thickness changes during film thinning. In contrast, for the native oxide, the position of the 
maximum stays at x = 6 mm, y = 4 mm until 187 h after preparation. By this time, the film is flat 
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Figure S2: Height profiles at constant j-position of TEHOS films shown in ?? during film thinning 
for films on (a) thermal (y = 5 mm) and (b) native oxide (y = 4 mm). 


o 

as can be seen from the film profiles shown in S2. With a thickness below 35 A, it is in the regime 
of molecular layering, as can bee seen from the thinning rates in ?? (a,b). For the flat film after 
200 h, the maximum values are found at the edge of the sample, where touching the sample with 
pincers for mounting and demounting plays a role. For the flat film after 211 h from preparation 
the /iniax are therefore shown in brackets. 

Table S2 : Maximum thickness (value and position) of TEHOS films shown in ?? during film 
thinning. Values in brackets are obtained at the edges within a flat film geometry. 


TEHOS film on thermal oxide TEHOS film on native oxide 


time [h] /imax [A] x [mm] y [mm] time [h] h^ax [A] x [mm] y [mm] 


20 

133.4 

3 

3 

18 

73.8 

6 

4 

44 

149.6 

6 

3 

43 

63.3 

6 

4 

69 

108.2 

6 

4 

66 

54.3 

6 

4 

91 

99.5 

6 

4 

89 

48.4 

6 

4 

116 

93.2 

5 

5 

114 

43.0 

6 

4 

140 

86.7 

5 

5 

138 

43.5 

6 

4 

164 

80.7 

5 

5 

162 

39.3 

6 

4 

188 

73.1 

5 

5 

187 

35.0 

6 

4 

212 

67.7 

5 

5 

211 

(33.9) 

9 

2 

236 

60.7 

5 

5 

235 

(34.7) 

9 

2 

260 

54.5 

4 

6 

259 

(34.8) 

9 

2 

284 

49.5 

4 

6 

283 

(31.1) 

1 

1 

308 

45.4 

6 

7 

306 

(32.1) 

9 

2 

332 

41.9 

6 

7 

331 

(28.1) 

1 

1 

356 

38.5 

5 

7 

354 

(22.3) 

9 

2 
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Thinning rates including subsets 


Average thinning rates calculated from different time spans At of the total time of investigation for 
the films on (a) native and (b) thermal oxide shown in ?? (a,c) native (A) and thermal (•). For 
native oxide (S3 a), for h < 40 A a scatter of the subsets about the thinning rate from the full data 
set is seen. This is due to large errors produced by averaging over not enough data by the particular 
binning method. Increased time of measurement produces more individual rates r, from thickness 
intervals h, 2 — (J/,! yielding enough data for the relevant bins of thickness h also for smaller h. 



Figure S3: Film thinning rates for TEHOS films on (a) native oxide (A) and (b) 100 nm thermal 
oxide (•) also shown in ?? (a,c); including thinning rates calculated from several subsets (differing 
in colors and symbols) of the full data sets. The subsets contain increasing numbers of succeeding 
measurements after sample preparation, and thus cover increasing time spans At of the total time 
of investigation. 

For the sample on thermal oxide (S3 b), the thickness dependent thinning rates decrease with 
increasing time span At in the thickness range between 80 and 40 A, as is also seen for the other 
film on thermal oxide shown in ?? (d). 


Derivation of interaction potentials from dielectric properties 

An approximation of van der Waals interactions may be obtained using the dielectric permittivities 
e, and refractive indices «, in the visible of the relevant materials (denoted by index i) together with 
the main electronic absorption frequency Vg in the UV, as long as similar Vg can be assumed for the 
used materials.^ The effective Hamaker coefficient A 12 - 23 , describing the strength of interaction 
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of material 1 with material 3 across medium 2 (for example Si02 with air across TEHOS) can be 
approximated by^ 


^ 12-23 


4 Vei+e2 


£ 3 -£ 2 \ 

£3+£2/ 


+ 



2 

2 


(S6) 


The first term in the sum gives the contribution from permanent dipoles, while the second term is 
related to dispersion. Effective interface potentials 4>(fi) for thin liquid films of thickness h can 
be modeled as sum of short-range steric repulsion C/h^ and the long-range non-retarded van der 
Waals potential ^vdw{h) according to^°’^^ 


^{h) 


C 


+ ^vdW (h) ■ 


(S7) 


In case of partially wetting polymer films, the strength C of the steric repulsion can be derived 
from fits to experimentally obtained effective interface potentials. 

Effective Hamaker coefficients for the used material compositions were calculated from optical 
properties using ??. S3 lists the obtained values, together with the corresponding contributions 
from permanent dipoles and from pure dispersion. Effective interface potentials for TEHOS 
films are calculated using ??. Thereby, the strength of the short-range repulsion C is estimated to 
be about 10“^^ Jm“^, since the interaction of the non-polar TEHOS with the polar Si02 is expected 
to be weaker than that of polar polystyrene with polar Si02 (10“^^ Jm“^), but stronger than the 


Table S3 : Effective Hamaker coefficients [zJ = 10 J] obtained using ??, and components from 
permanent dipoles and from dispersion. 


material composition 

permanent dipole part 

dispersion part 

effective Hamaker coefficient 

Si/TEHOS/air 

-0.7 

-41.7 

-42.5 

Si02/TEHOS/air 

0.2 

-2.4 

-2.3 

Si/hexane/air 

-0.7 

-37.1 

-37.7 

Si02/hexane/air 

0.1 

-9.0 

-8.8 

Si/TEHOS/hexane 

-0.1 

-6.1 

-6.2 

Si02/TEHOS/hexane 

0.02 

-0.3 

-0.3 
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interaction of polystyrene with the non-polar OTS (10“^^ For the silicon substrate with 

100 nm thermal oxide, the simplified ?? can be used with Asi02 to calculate ^^dwih), because 
for a thickness ^/ = 100 nm of the coating oxide layer, the contribution from Asi is negligible in 
??. The average thickness of native oxide was determined by ellipsometry to be 2 nm. In case of 
native oxide, ?? has to be used with a thickness d = 2 nm of the coating layer. This yields a notable 
contribution from the underlying silicon. 
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